home *** CD-ROM | disk | FTP | other *** search
/ SPACE 2 / SPACE - Library 2 - Volume 1.iso / program / 316 / libsrc / fixnum.s < prev    next >
Encoding:
Text File  |  1988-10-20  |  4.7 KB  |  250 lines
  1. # /* Subroutines needed by GCC output code on some machines.  */
  2. # /* Compile this file with the Unix C compiler!  */
  3. # 
  4. # #include "config.h"
  5. # 
  6. # union double_di { double d; int i[2]; };
  7. # union flt_or_int { int i; float f; };
  8. # 
  9. # #ifdef WORDS_BIG_ENDIAN
  10. # #define HIGH 0
  11. # #define LOW 1
  12. # #else
  13. # #define HIGH 1
  14. # #define LOW 0
  15. # #endif
  16. # 
  17. # #ifdef L_umulsi3
  18. # _umulsi3 (a, b)
  19. #      unsigned a, b;
  20. # {
  21. #   return a * b;
  22. # }
  23. # #endif
  24.  
  25. .text
  26.     .even
  27. .globl __umulsi3
  28. __umulsi3:
  29.     link a6,#0        | 0 bytes of local
  30.     movel d2,sp@-        | push d2
  31.     movel a6@(8),d2        | get a
  32.     beq umul_zero        | 0? return 0
  33.     movel a6@(12),d1    | get b
  34.     beq umul_zero        | 0? return 0
  35.  
  36.     mulu d1,d2        | mul bottom parts
  37.     movel d2,d0        | save that
  38.     movew a6@(8),d2        | get hi part of a
  39.     beq umul_1        | zero, skip it
  40.     movel a6@(12),d1    | get b
  41.     mulu d1,d2        | mul a hi by b lo
  42.     lsll #8,d2        | shift over
  43.     lsll #8,d2        |  ...twice...
  44.     addl d2,d0        | add that in
  45. umul_1:
  46.     clrl d1
  47.     movew a6@(12),d1    | get b hi
  48.     beq umul_ret
  49.     movel a6@(8),d2        | and a lo
  50.     mulu d1,d2        | mult those
  51.     lsll #8,d2        | shift it over
  52.     lsll #8,d2        |  ... twice...
  53.     addl d2,d0        | and add that in
  54.     bra umul_ret        | go home
  55. umul_zero:
  56.     clrl d0            | return 0
  57. umul_ret:
  58.     movel sp@+,d2        | get d2 back
  59.  
  60.     unlk a6            | flush frame
  61.     rts
  62. #
  63. # #ifdef L_mulsi3
  64. # _mulsi3 (a, b)
  65. #      int a, b;
  66. # {
  67. #   return a * b;
  68. # }
  69. # #endif
  70.  
  71. .text
  72.     .even
  73. .globl __mulsi3
  74. __mulsi3:
  75.     link a6,#-2        | one word of local
  76.     clrw a6@(-2)        | zap the local; it's our neg flg
  77.  
  78.     movel a6@(12),d0    | get b
  79.     bpl muls1        | pos, ok
  80.     negl d0            | negate it
  81.     addqw #1,a6@(-2)    | and inc the flag
  82. muls1:
  83.     movel d0,sp@-        | push adjusted b
  84.     movel a6@(8),d0        | get a
  85.     bpl muls2        | pos, ok
  86.     negl d0            | negate it
  87.     subqw #1,a6@(-2)    | and dec the flag
  88. muls2:
  89.     movel d0,sp@-        | push adjusted a
  90.     jbsr __umulsi3        | do an unsigned mult
  91.     addql #8,sp        | flush args
  92.     tstw a6@(-2)        | flag set?
  93.     beq muls3        | nope, go ahead and return
  94.     negl d0            | negate the result
  95. muls3:
  96.  
  97.     unlk a6            | unlink the frame
  98.     rts            | done!
  99. #
  100. # #ifdef L_udivsi3
  101. # _udivsi3 (a, b)
  102. #      unsigned a, b;
  103. # {
  104. #   return a / b;
  105. # }
  106. # #endif
  107.  
  108. .text
  109.     .even
  110. .globl __udivsi3
  111. __udivsi3:
  112.     link a6,#0
  113.  
  114.     movel d2,sp@-        | check our bags
  115.     movel d3,sp@-
  116.     movel d4,sp@-
  117.  
  118.     clrl d0            | d0 will be accum
  119.     movel a6@(8),d1        | d1 is a, dividend
  120.     movel a6@(12),d2    | d2 is b, divisor
  121.     movel #1,d3        | get a 1 into bit 0 of bit mask
  122.     movel #31,d4        | d4 is bit num for testing need to sub
  123. udiv_sl:
  124.  
  125.     btst d4,d2        | bit get to top yet?
  126.     bne udiv_sr        | yes, start testing for subtraction
  127.     lsll #1,d3        | shift bitnum left
  128.     lsll #1,d2        | shift divisor left
  129.     bra udiv_sl
  130. udiv_sr:
  131.     cmpl d2,d1        | divisor > dividend?
  132.     beq udiv_sub
  133.     bcs udiv_s
  134. udiv_sub:
  135.     subl d2,d1        | take divisor (shifted) out of dividend
  136.     orl d3,d0        | set this bit in d0
  137.     tstl d1            | anything left?
  138.     beq udiv_done        | nope
  139. udiv_s:
  140.     subql #1,d4        | dec bitnum (new only)
  141.     lsrl #1,d2        | shift right one
  142.     lsrl #1,d3        | and the bitnum too
  143.     bne udiv_sr        | more bits, keep dividing
  144.  
  145. udiv_done:
  146.     movel sp@+,d4
  147.     movel sp@+,d3        | get bags back
  148.     movel sp@+,d2
  149.  
  150.     unlk a6
  151.     rts
  152. #
  153. # #ifdef L_divsi3
  154. # _divsi3 (a, b)
  155. #      int a, b;
  156. # {
  157. #   return a / b;
  158. # }
  159. # #endif
  160.  
  161. .text
  162.     .even
  163. .globl __divsi3
  164. __divsi3:
  165.     link a6,#-2        | 2 bytes of local
  166.     clrw a6@(-2)        | zap the local; it's our neg flg
  167.  
  168.     movel a6@(12),d0    | get b
  169.     bge divs1        | pos, ok
  170.     negl d0            | negate it
  171.     addqw #1,a6@(-2)    | and inc the flag
  172. divs1:
  173.     movel d0,sp@-        | push adjusted b
  174.     movel a6@(8),d0        | get a
  175.     bge divs2        | pos, ok
  176.     negl d0            | negate it
  177.     subqw #1,a6@(-2)    | and dec the flag
  178. divs2:
  179.     movel d0,sp@-        | push adjusted a
  180.     jbsr __udivsi3        | do an unsigned div
  181.     addql #8,sp        | flush args
  182.     tstw a6@(-2)        | flag set?
  183.     beq divs3        | nope, go ahead and return
  184.     negl d0            | negate the result
  185. divs3:
  186.  
  187.     unlk a6            | unlink the frame
  188.     rts            | done!
  189. #
  190. # #ifdef L_umodsi3
  191. # _umodsi3 (a, b)
  192. #      unsigned a, b;
  193. # {
  194. #   return a % b;
  195. # }
  196. # #endif
  197.  
  198. .text
  199.     .even
  200. .globl __umodsi3
  201. __umodsi3:
  202.     link a6,#-4
  203.     movl d2,sp@-
  204.     movel a6@(12),sp@-
  205.     movel a6@(8),sp@-
  206.     jbsr __udivsi3
  207.     movel d0,d0
  208.     movel a6@(12),sp@-
  209.     movel d0,sp@-
  210.     jbsr __mulsi3
  211.     movel d0,d2
  212.     movel a6@(8),d1
  213.     subl d2,d1
  214.     movel d1,a6@(-4)
  215.     movel a6@(-4),d0
  216.     movl a6@(-8),d2
  217.     unlk a6
  218.     rts
  219. # 
  220. # #ifdef L_modsi3
  221. # _modsi3 (a, b)
  222. #      int a, b;
  223. # {
  224. #   return a % b;
  225. # }
  226. # #endif
  227.  
  228. .text
  229.     .even
  230. .globl __modsi3
  231. __modsi3:
  232.     link a6,#-4
  233.     movl d2,sp@-
  234.     movel a6@(12),sp@-
  235.     movel a6@(8),sp@-
  236.     jbsr __divsi3
  237.     movel d0,d0
  238.     movel a6@(12),sp@-
  239.     movel d0,sp@-
  240.     jbsr __mulsi3
  241.     movel d0,d2
  242.     movel a6@(8),d1
  243.     subl d2,d1
  244.     movel d1,a6@(-4)
  245.     movel a6@(-4),d0
  246.     movl a6@(-8),d2
  247.     unlk a6
  248.     rts
  249.  
  250.